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Abstract 

Background: Periodontitis is an infectious and inflammatory disease of polymicrobial etiology that can lead to the 
destruction of bones and tissues that support the teeth. The management of chronic periodontitis (CP) relies 
heavily on elimination or at least control of known pathogenic consortia associated with the disease. Until now, 
microbial plaque obtained from the subgingival (SubG) sites has been the primary focus for bacterial community 
analysis using deep sequencing. In addition to the use of SubG plaque, here, we investigated whether plaque 
obtained from supragingival (SupG) and tongue dorsum sites can serve as alternatives for monitoring CP-associated 
bacterial biomarkers. 

Results: Using SubG, SupG, and tongue plaque DNA from 1 1 healthy and 13 diseased subjects, we sequenced 
V3 regions (approximately 200 bases) of the 16S rRNA gene using lllumina sequencing. After quality filtering, 
approximately 4.1 million sequences were collapsed into operational taxonomic units (OTUs; sequence identity 
cutoff of >97%) that were classified to a total of 19 phyla spanning 1 14 genera. Bacterial community diversity and 
overall composition was not affected by health or disease, and multiresponse permutation procedure (MRPP) on 
Bray-Curtis distance measures only supported weakly distinct bacterial communities in SubG and tongue plaque 
depending on health or disease status (P<0.05). Nonetheless, in SubG and tongue sites, the relative abundance of 
Firmicutes was increased significantly from health to disease and members of Synergistetes were found in higher 
abundance across all sites in disease. Taxa indicative of CP were identified in all three locations (for example, 
Treponema denticola, Porphyromonas gingivalis, Synergistes oral taxa 362 and 363). 

Conclusions: For the first time, this study demonstrates that SupG and tongue dorsum plaque can serve as 
alternative sources for detecting and enumerating known and novel bacterial biomarkers of CP. This finding is 
clinically important because, in contrast with SubG sampling that requires trained professionals, obtaining plaque 
from SupG and tongue sites is convenient and minimally-invasive and offers a novel means to track CP-biomarker 
organisms during treatment outcome monitoring. 
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Background 

Bacterial communities that colonize the human body are 
intimately linked with host physiology, immunity, me- 
tabolism, and nutrition [1-3]. Of colonization sites, the 
mouth accommodates a bacterial consortium with high 
taxonomic richness and phylogenetic diversity [4-6]. The 
microbial composition of the plaque biofilm has a critical 
role in oral health. Disruption of plaque homeostasis can 
stimulate tissue destruction and inflammation, leading to 
infections such as dental caries, gingivitis, and periodon- 
titis [4,7]. A comprehensive and systematic profiling of the 
oral biofilm is necessary, not only to understand microbial 
associations with localized infections, but also because the 
oral microbiome has long been known as a reservoir for 
infections at other body sites [8,9]. Bacterial species nor- 
mally found in the oral cavity have been associated with 
distal infections in the lungs, heart, brain, and liver, either 
reflecting their involvement in opportunistic infections 
as a result of systemic changes in the body, or perhaps 
suggesting a causative link [8,10-13]. Therefore, the plaque 
consortia can be a mirror and possibly even a monitor of 
oral and non-oral health and disease. 

Periodontitis is a chronic inflammatory disease affecting 
tooth-supporting structures including the alveolar bone, 
connective tissue attachment, and gingiva [14,15]. Although 
several forms of periodontal diseases have been recognized, 
the predominant category is chronic periodontitis (CP), 
which remains a primary cause of tooth loss in adults 
worldwide [16]. In addition to its large socio-economic 
burden [17], CP is associated with considerable morbidity 
in terms of pain, uncomfortable chewing, oral malodor, and 
tooth migration. Although the pathogenesis of periodontitis 
is multifactorial and includes genetic and epigenetic factors, 
the development of periodontitis is modulated by microbial 
biofilm that forms on and around teeth, eliciting an inflam- 
matory host reaction [18-21]. The etiology of periodontitis 
is polymicrobial. In particular, disease progression has been 
linked with the proliferation of Gram-negative anaerobic 
species such as Porphyromonas gingivalis, Treponema 
denticola, and Tannerella forsythia [22]. 

Over the past two decades, culture-based, immuno- 
histochemistry, and molecular techniques have focused 
largely on a small subset of CP-associated microorgan- 
isms in the gingival sulcus and thus have overlooked 
the impact of potentially large numbers of other oral 
bacterial species in the infectious process [22-24]. Al- 
though these as-yet-unidentified organisms might not 
always be associated intimately with the periodontium 
at all times, it is still likely that several of them, at one 
time or another, participate in the pathophysiological 
processes that lead to and cause the advancement of 
CP. Emerging evidence of subgingival (SubG) community 
analysis using DNA sequencing suggests that periodontal 
destruction is associated with many uncultivable and 



uncharacterized disease indicator organisms [25-27]. How- 
ever, despite the growing list of periodontitis -associated 
bacteria, there is slow progress in the use of this informa- 
tion to develop more cost-effective and sensitive methods 
to treat, control, or prevent disease progression. It is prob- 
able that a more complete characterization of bacterial 
biomarkers of CP will lead to the development of new ther- 
apeutics, improved diagnostics, and alternative methods for 
monitoring the outcomes of treatment success. 

Recently, high-throughput sequencing approaches have 
assessed bacterial composition in deep and shallow SubG 
samples derived from patients with periodontitis compared 
to healthy patients [25]. The results demonstrated that pa- 
tients with periodontitis had more diverse combinations of 
species but also had all of the health-associated species, al- 
beit at a lower frequency. These results agreed with another 
study that showed increasing diversity in diseased compared 
to healthy sites [27]. However, diversity did not increase be- 
tween bleeding and non-bleeding SubG sites; an increase in 
total bacterial load was demonstrated in the SubG samples 
obtained from bleeding sites. Both of these studies noted a 
clear distinction between healthy and diseased SubG sites 
when using distance metric analysis [25,27]. 

In this study, we used 16S rRNA gene analyses for exam- 
ining bacterial communities associated with periodontal 
health and CP, using plaque extracted from SubG and 
supragingival (SupG) as well as sites on the dorsum of the 
tongue. We hypothesized that consistent differences in bac- 
terial communities of both dental and tongue plaque would 
exist between CP and healthy patients and that these differ- 
ences would be largely due to known disease-associated taxa 
as well as other as-yet-uncharacterized bacteria. Another 
goal was to determine if SupG and tongue plaque sampling 
could be used as a novel and non-invasive detection method 
of indicator taxa for CP, instead of conventionally sam- 
pled SubG plaque. A relatively simpler alternative sam- 
pling protocol would greatly facilitate the detection of 
periodontitis-associated bacterial biomarkers in disease 
monitoring. For the first time, we provide evidence that 
both SupG and tongue dorsum sites may be suitable as al- 
ternative sampling sources for the detection and enumer- 
ation of selected CP-associated bacterial biomarkers. The 
use of tongue plaque or SupG plaque that is above the gum 
line to monitor oral health is convenient and circumvents 
the utility of specially trained personnel for the sampling 
procedure; in fact, samples can be provided by the patients 
themselves. From a clinical standpoint, this study provides 
important groundwork for developing a minimally invasive 
method to monitor CP treatment outcomes. 

Results 

Demographics, clinical attributes, and sequencing 

Fifteen women and nine men in the age range of 18 to 
71 years were recruited for this study. The cohort with a 
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healthy periodontium consisted of eight women and 
three men in the age range of 18 to 56 years; the CP 
group comprised seven women and six men in the age 
range of 25 to 71 years (Table 1). CP was diagnosed on the 
basis of radiographic and clinical examination. Parameters 
derived from clinical evaluation included probing depth 
and clinical attachment loss (CAL). Bone loss on each tooth 
was assessed radiographically. Patients in the CP group 
had >30% sites with >5 mm of CAL and probing pocket 
depths of >6 mm, whereas the healthy group had >30% 
sites with <2 mm attachment loss and no pockets >3 mm. 
No patients with systemic diseases, antibiotic use, and den- 
tal cleanings within a 6-month time frame were included in 
the study. Differences in patient demographics and clinical 
data (Table 1) were assessed with a Welch's two sample 
t test or a chi-squared test of independence. No significant 
differences in age, weight, gender, smoker status, or drinker 
status were found between health and CP. All patient and 
sample metadata are available in Additional file 1. 

For comparing bacterial community composition be- 
tween participants with a healthy periodontium and 
CP patients, DNA isolated from plaque was subjected 
to PCR amplification using primers covering an ap- 
proximate 200-bp segment of the V3 region of the 16S 
rRNA gene. DNA sequencing yielded approximately 
5.8 million sequences, with a minimum of 27,187 se- 
quences per sample (average of 57,441 ± 23,336) and 
between-sample analyses rarefied all sample data to 
27,187 sequences. After quality filtering and PANDAseq 
assembly (minimum overlap of 60 nucleotides), approxi- 
mately 4.1 million sequences remained. These sequences 
were collapsed into operational taxonomic units (OTUs) 
based on a pairwise sequence identity cutoff of 97%. Of the 
total sequences, 49,327 formed singletons and 6,160 OTUs 
did not classify to the phylum level (representing 66,466 



Table 1 Demographics and clinical parameters of participants 



Demographics, clinical 
parameters, and habits 


Healthy group 
(n = 11) 


Diseased group 
(n = 13) 


Sex 


3 men; 8 women 


6 men; 7 women 


Age (years) 


38.4 ±4.1 


46.8 ± 4.0 


Weight (lbs) 


146.9 ±7.1 


159.1 ±8.2 


Probing depth (mm) 
in > 30% of sites 


<2 


>5 


Attachment loss (mm) 
in >30% of sites 


<1 


>3 


Periodontal diagnosis 3 


Periodontal health 


Chronic periodontal 
disease (CP) 


Smokers b 


1 


6 


Drinkers c 


2 


5 



Periodontal diagnosis based on criteria as outlined in the 1999 American 
Academy of Periodontology Diagnostic Guidelines; b a smoker was defined as 
someone who had one or more cigarettes per day for >3 months; c a drinker 
was defined as an individual who consumed >56 drinks per year. 



sequences, approximately 0.02% of the usable sequences). 
OTUs were classified to a total of 19 phyla, 114 genera, 
and 302 species. 

Advanced periodontal disease is not correlated with 
changes in a-diversity 

Recent studies aimed at assessing site-specific bacterial 
community diversity have implicated considerable changes 
in the taxonomic composition in health versus disease 
[25,27,28]. Hence, we examined whether the biofilm 
community composition in the CP cohort could be 
correlated with consistent alterations in OTU diversity. 
The a-diversity within SupG, SubG, and tongue sample 
sites was investigated using Phyloseq. Standard diversity 
metrics were evaluated, including observed richness, 
Shannon index, and Simpson index. There were no sta- 
tistically significant differences in a-diversity between 
healthy and diseased individuals across all sampling 
sites and metrics (Figure 1). The average a-diversity scores 
for all metrics were higher in disease for SupG and SubG 
sites, but lower in disease for the tongue. 

Tongue harbors a unique consortium, and periodontal 
destruction is associated with a weakly distinct shift 
relative to health 

To assess structural similarities in the biofilm communities 
among SupG, SubG, and tongue sample sites in healthy 
and diseased patients, a multiresponse permutation proced- 
ure (MRPP) analysis based on a Bray-Curtis distance matrix 
was used. To visualize the results, a non-metric multidi- 
mensional scaling (NMDS) plot (Figure 2) and a principal 
coordinates analysis (PCoA) plot were generated (Figure 3), 
both based on Bray-Curtis distances. Bray-Curtis provides a 
measure of community composition differences between 
samples based on OTU counts, regardless of taxonomic 
assignment. Ordinations based on this metric demon- 
strated a clear separation of tongue samples from tooth 
sample locations, indicating that the bacterial community 
on the tongue is distinct from those in SupG and SubG 
sites (A = 0.06, T= -22.7, P< 0.001). Alternatively, there 
was no obvious separation between SupG and SubG sam- 
ples. Similar results were shown using a PCoA plot based 
on UniFrac distances (plot available in Additional file 2). 

Further investigation was carried out to determine 
whether or not there are distinct bacterial compositions 
in SupG, SubG, and tongue sites depending on the 
health status of the periodontium. For each site, a PCoA 
ordination and MRPP analysis, based on Bray-Curtis 
distances, was performed. The PCoA ordination did not 
reveal strong grouping of health and disease on the primary 
axes (Figure 3). In addition, no clear separation of health 
and disease samples was seen in the UniFrac-based PCoA 
plots (see Additional file 2). MRPP analysis of the distance 
matrix suggested very weak differences between health and 
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Figure 1 Alpha diversity calculations for healthy (red) and diseased (blue) samples. No significant differences between healthy and 
diseased samples were observed for various a-diversity measures. Calculations and plots were generated using the Phyloseq software package. 
(A) All samples (SubG, SupG, and Tongue) pooled. (B) SubG samples. (C) SupG samples. (D) Tongue samples. 




Figure 2 Bray-Curtis based non-metric multidimensional scaling 
(NMDS) plot of all samples. NMDS plot (stress 0.22) shows the 
clear clustering of tongue samples distinct from SubG and SupG 
samples. Ellipsoids represent a 95% confidence interval surrounding 
each group. MRPP analysis concluded that the members of the 
tongue and SubG/SupG groups were more dissimilar than expected 
by chance {A = 0.06, T= -22.7, P < 0.001). 



disease samples for the combined, SubG, and tongue sites 
CP < 0.05; full MRPP results available in Additional file 3), 
although not for the SupG site. 

Relative abundance of distinct taxa is affected by 
periodontal destruction in oral plaque communities 

Using CP and healthy samples, the relative abundance of 
taxa were compared to determine whether periodontal 
destruction was correlated with substantial changes in the 
abundance of specific bacterial taxa. The abundances of 
each taxonomic group were assessed for significant differ- 
ences between health and disease with a Mann- Whitney 
test. The most abundant taxa in CP-associated biofilms 
belonged to the phyla Bacteroidetes and Proteobacteria; 
the most significant difference occurred in the phylum 
Synergistetes (Figure 4). To define an association with 
SupG, SubG, and tongue locations, the relative abundance 
of taxa by sample site was also examined. In SubG and 
tongue samples, the relative abundance of Firmicutes was 
significantly increased under CP versus health (Figures 4B, 
D). Despite being associated with CP subjects in the sample 
site combined analysis, Bacteroidetes and Proteobacteria 
were found to be associated with healthy subjects in SubG 
samples (Figure 4B). Across all sites, the Synergistetes 
phylum was found in higher abundance in patients with 
CP as compared to healthy subjects (Figure 4). 
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Figure 3 Principal coordinate analysis (PCoA) plot with Bray-Curtis dissimilarity. Results revealed that tongue samples clustered separately 
from tooth samples (SubG and SupG), suggesting that the tongue community is relatively unique from that in SupG and SubG sites. No obvious 
clustering was apparent within or between SupG and SubG samples. (A) Bray-Curtis PCoA ordination of pooled SupG, SubG, and tongue samples. 
(B) Ordination of SubG samples. Healthy patient samples are shown in green, diseased patient samples in pale red. (C) Ordination of SupG 
samples. (D) Ordination of tongue samples. 



Closer examination of CP-associated abundance profiles 
at a higher taxonomic resolution revealed many highly 
abundant genera in the SubG to be significantly increased 
in health, including Corynebacteriurn, Capnocytophaga, 
Campylobacter, Neisseria, and Kingella (shown in 
Additional file 4, panel A). On the other hand, the 
tongue harbored many highly abundant disease-associated 
genera such as Treponema, Synergistes, and Clostridials 
(shown in Additional file 4, panel C). Examining results at 
the species level highlighted classic periodontal disease- 
associated OTUs in each of the three sites. These in- 
cluded T. forsythia, Filifactor alocis, and Porphyromonas 
endodontalis (see Additional file 5). In particular, F. alocis 
was found in all three sites in high abundance and signifi- 
cantly higher proportions in disease. Health-associated taxa 
included unclassified Fusobacteria and Corynebacteriurn 
matruchotii in the SubG site, Granulicatella adiacens 
and Selenomonas artemidis in SupG, and S. artemidis 
and unclassified Porphyromonas in tongue samples 
(significantly elevated species listed in Additional file 6). 



Periodontal destruction and periodontal health are 
associated with indicator organisms 

To determine if organisms present in the sample sites can 
serve as specific indicators of health or disease, indicator 
analysis of the OTUs was performed. OTUs were scored by 
their abundance and presence/absence in the health and 
disease groups. Indicator OTUs were those which were sig- 
nificantly more abundant and present in all samples be- 
longing to one group and also absent or low abundance in 
the other group [29]. At the species level, SubG, SupG, and 
tongue sites of the CP cohort were associated with 10, 11, 
and seven species, respectively, whereas 13 and six were 
health-associated indicators for SubG and SupG sites, re- 
spectively (Figure 5). No indicators of health were found on 
the tongue at the specified indicator value threshold and 
abundance. In all sites, the vast majority of OTUs were not 
significant indicators of either health or disease. 

Indicator species analysis of combined data for SupG, 
SubG, and tongue sites showed several indicator OTUs as 
strong representatives of either disease or health states 



Galimanas et al. Microbiome 2014, 2:32 
http://www.microbiomejournal.eom/content/2/1/32 



Page 6 of 13 



0.100 ■ 



0.001 



** 



ffl 



• i 



I I 



W nP" 

J> a? & 



4" 



& c$ 



0.100 - 



0.001 - 









** 

- =^=| ** 

m X . 


B 




! 


^ -r- 










I sk* 











1e-02 



1e-03 ■ 



1e-04 ■ 



D 



0.100 - 



0.001 - 




Phylum 

Figure 4 Health and disease associated phyla, by relative 
abundance. The top five most abundant phyla with significant 
differences between health and disease are shown. Significance was 
assessed with a Mann-Whitney test (P< 0.05). Disease abundances 
are shown in pale red, and health abundances in green. (A) All samples 
(SubG, SupG, and tongue). (B) SubG samples. (C) SupG samples. 
(D) Tongue samples. *0.01 < P < 0.05; **0.001 < P < 0.01 ; 
***P < 0.001 ; ****p< 0.0001. 



(Figure 6A). Synergistes oral taxon (OT) 363 and Synergistes 
OT 362 showed the highest indicator values for CP; 
Streptococcus sanguinis, T. forsythia, F. alocis, Dialister 
invisus, Streptococcus sp., and TM7 401H12 were also 
strong indicators for CP (see Additional file 6). F. alocis was 
a strong indicator of disease in SubG and SupG samples, 
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Figure 5 Numbers of taxa associated with health and disease 
by indicator species analysis. Dufrene-Legendre indicator species 
analysis revealed OTU associations with either health or disease. 
The number of indicator taxa at phylum, genus, and species levels 
(OTU indicator value >0.5 and taxa sequence abundance >1 00) of 
health (green), disease (pale red), or both (light blue) are shown. 
Taxa that were present in the samples but not an indicator of either 
health or disease are shown in grey. (A) SubG samples. (B) SupG 
samples. (C) Tongue samples. 



but a relatively weak indicator of disease on the tongue 
(indicator value of 0.46; not shown in Figure 7). Interest- 
ingly, T. forsythia, as well as the unclassified Synergistes 
OT 363 and 362, were indicators of disease in this in- 
vestigation (Figure 6). Health indicators included un- 
classified Fusobacteriales, TM7 7BB428, Campylobacter 
rectus, uncultured Lachnospiraceae, OT 100, unclassified 
Bacteroidetes, Aggregatibacter, Capnocytophaga gingivalis, 
C. matruchotii, Neisseria, Selenomonas noxia, and 
Selenomonas artemidis. A full list of indicator taxa and 
related statistics is available in Additional file 7. 

Quantitative PCR confirms periodontal disease indicators 

There were no significant alterations in bacterial 16S rRNA 
gene copy numbers (per ng of DNA) in healthy versus 
diseased plaque (Figure 7). Using species-specific primers, 
qPCR was used to determine the relative abundance of 
three species (that is, P. gingivalis, F. alocis, and T. forsythia) 
in SubG, SupG, and tongue sites using healthy and diseased 
plaque. Of these species, P. gingivalis and T. forsythia were 
validated as the most significantly abundant taxa in the CP 
group (Figure 7), in agreement with the sequence data. In 
particular, levels of P. gingivalis were elevated in all three 
sites affected by CP, suggesting that SupG and tongue are 
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Figure 6 Bubble plot of indicator OTUs associated with health and disease. Each bubble represents an OTU identified by Dufrene-Legendre 
indicator species analysis as being associated with health (green) or disease (pale red). (A) Pooled SupG, SubG, and tongue samples; (B) SubG 
samples; (C) SupG samples; and (D) tongue samples. Bubble area is proportional to OTU sequence abundance. Only OTUs with indicator values >0.5 
and taxa with a total of >1 00 sequences are shown. Complete indicator OTU lists are available in Additional file 7. 



as effective as SubG sampling sites for quantification of 
P. gingivalis (Figure 7). The relative number of T.forsythia 
was significantly elevated in diseased samples in SupG and 
SubG sites, but not in the tongue. The relative abundance 
of F. alocis was elevated significantly in SupG diseased 
samples, relative to controls. Although its abundance was 
markedly increased in diseased SubG samples, this result 
could not be validated on the basis of statistical analyses 
performed in this study. Some disagreement between 
the high-throughput sequencing and qPCR was seen. 
The sequencing results showed a significant increase 
in F. alocis in the disease state in all three sites, but the 
qPCR results suggested an increase in F. alocis in the 
healthy samples of the tongue site (qPCR data shown 
in Figure 7, Additional file 5 visualizes the differences 
between species, and Additional file 6 lists species with 
significantly different abundances and relevant statistics). 
An increase in P. gingivalis was indicated in all three 
sites by the qPCR method, but this increase could not 



be validated with statistical analyses of the sequence 
data. A significant increase in T.forsythia was reported 
in the SubG and SupG sites by qPCR, but sequencing 
results identified an increase in this organism in the 
SubG and tongue sites. 

Discussion 

In this study, microbial communities associated with CP 
and periodontal health were characterized by examining 
plaque bacteria from SupG, SubG, and tongue sites, 
providing a comprehensive cross-sectional description 
of microbiota in the presence of CP and periodontal 
health. Findings reported here show that most OTUs are 
shared between periodontal health and disease. Nonetheless, 
a relatively small proportion of OTUs were distinct to dis- 
ease and included T. denticola, T. forsythia, P. endodontalis, 
Synergistes OT 362, and Synergistes OT 363. Because these 
disease indicators were present in high abundance in tongue 
plaque of individuals suffering from CP, we propose that 
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within each of the respective locations in the oral cavity is based on non-parametric Mann-Whitney Test with outliers removed using Grubbs' 
Outlier Test. *0.01 < P < 0.05; **0.001 < P < 0.01 ; ***P < 0.001 ; ****P < 0.0001 . 



tongue plaque may be used as an easily accessible, non- 
invasive source to collect plaque to quantify bacterial 
biomarkers to monitor treatment outcomes against CP 
over time. In developing such methodologies, it will be 
useful to evaluate the ratio of selected periodontal 
disease-associated taxa to health- associated taxa and with 
larger numbers of human participants for validation. It 
will also be critically important to differentiate between 
active and inactive periodontal disease because the micro- 
bial communities under these two conditions can vary. 

High-throughput sequencing studies conducted to date 
have examined plaque derived from SubG and SupG sites 
to differentiate between healthy and periodontal disease 
consortia [25,27,30]. SubG samples can be more difficult 
to obtain as compared to samples taken from SupG or 
tongue sites. Collecting plaque from SubG locations 
generally requires that dentally trained investigators 
do most of the sampling. However, if sampling of oral 



microbes associated with periodontitis were to be 
done on a larger basis and for large population sam- 
ples (for example, investigations that seek to correlate 
periodontitis with non-oral diseases), it would be most 
advantageous to develop a simpler approach to microbial 
sampling. Here, we investigated the use of SupG and 
tongue samples as potential surrogate sites for bacterial 
biomarker detection under health and CP. The keratinized 
mucosal surface present in the tongue dorsum is a distinct 
colonization surface relative to that of tooth structures, 
which are comprised of highly mineralized tissues such as 
enamel and cementum. Tooth surfaces are colonized by 
different consortia depending on the distinct anatomical 
site and gingival surfaces [31,32]. Not surprisingly, in our 
ordination plots using SupG, SubG, and tongue plaque, 
tongue samples clearly separated from tooth plaque 
samples (Figures 2 and 3A). This separation indicated 
that bacterial communities present in tongue sites were 
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distinct from those found in SupG and SubG sites. 
Previously, Griffen and colleagues reported clustering 
based on disease, although only SubG plaque was sampled 
in their study [25]. Similar clustering was not seen in our 
PCoA ordinations (Figure 3). Similar to the weak separation 
found in our findings, a recent metagenomic and metatran- 
scriptomic study of SubG and SupG plaque failed to find a 
strong separation between health and CP [33]. The shift 
from health to disease may be due to a subtle ecological 
shift and as a result will not appear in the primary axes of 
the ordination plots. An MRPP analysis, which calculates 
group separation based on higher dimensional data, would 
be able to identify a more subtle separation by utilizing 
all dimensions of the data. MRPP revealed small but 
significant differences in the community composition 
between health and diseased states in the combined, 
SubG, and tongue sites (P < 0.05; full MRPP statistics 
available in Additional file 3). 

Previous work demonstrated that certain bacterial com- 
plexes associate with each other and potential periodontal 
diseases [22]. DNA-DNA hybridization was used to moni- 
tor 40 bacterial species using 13,261 plaque samples from 
185 healthy or periodontal disease patients. By examining 
pre- and post-therapy bacterial profiles, species prevalence, 
community ordination, and association with 0 to 3 mm, 4 
to 6 mm, and >6 mm probing depths, disease-associated 
bacteria included P. gingivalis, T.forsythia, and T. denticola, 
which were highly prevalent in deeper sites that bled upon 
probing [34]. Prior to this, it was reported that that the 
probability of finding oral treponemes at sites that harbored 
P. gingivalis was escalated as periodontal destruction was 
increased [35]. Others demonstrated that species present in 
disease were also present in the healthy state, albeit in lower 
proportions [25], while another group demonstrated that 
SupG sites harbored the same species as SubG sites below 
the gingival margin, again in smaller proportions [36]. The 
latter study postulated that this may be conducive to re- 
infection and individuals with 'healthy sites with abundant 
disease-indicator organisms may be at higher risk for dis- 
ease progression. As a result, it may not be unexpected that 
the Bray-Curtis and UniFrac distance ordinations did not 
show consistent group separation. Furthermore, additional 
statistical tests of the datasets may be more informative 
than ordinations for identifying disease-associated taxa. 

In addition to classic periodontitis-associated pathogens 
mentioned above, we and others detected F. alocis in 
patients with increased periodontal destruction (Figure 6) 
[37-39]. In this study, F. alocis was detected in signifi- 
cantly higher abundances in disease in all three sites 
(differences between health and disease are visualized in 
Additional file 5, and statistics are available in Additional 
file 6). F. alocis is a rod-shaped Gram-positive facultative 
anaerobe that possesses trypsin proteases, invades human 
cells, resists oxidative stress, and forms biofilms, which 



are properties conducive to disease initiation and progres- 
sion [38,40,41]. In particular, the prevalence of F alocis was 
elevated in patients suffering from generalized aggressive 
periodontitis and CP, whereas it was rarely detected in the 
control group resistant to periodontitis [37]. Although this 
pathogen was shown to colonize apical parts of the pocket 
in close proximity to the soft tissues, it was also shown 
to contribute to the structural integrity of these multi- 
species biofilms. Hence, F. alocis was suggested to be a 
good diagnostic marker organism for periodontal dis- 
ease [37]. Based on previous investigations, as well as 
the data shown here, here, we propose that F. alocis be 
included as a novel member in Socranskys 'Red complex' 
of periodontal pathogens [22], 

Comparisons within sites using the qPCR assay demon- 
strated an increase in bacterial abundance only in the 
tongue samples derived from patients with CP (Figure 7). 
The total bacterial abundance was not significantly altered 
when SubG and SupG sites were compared. This result 
is in contrast to recent studies that indicated a marked 
increase in the total bacterial abundance in SubG sites of 
subjects affected by periodontal disease [36,42]. Our qPCR 
analysis, as well as those conducted by others, does not 
account for the total bacterial biomass of plaque samples. 
As in other studies, the total bacterial load between health 
and disease states were calculated by normalizing qPCR 
results with the total nanogram of genomic DNA in each 
sample. Hence, it is not possible to make inferences about 
the total bacterial load between health and disease sam- 
ples (that is, per mg plaque dry weight) using the methods 
employed in this study. Further, although qPCR was useful 
for validating our Illumina-based indicator species and 
results generally agreed, observed trends could not al- 
ways be statistically validated by both methods. In the 
case of F. alocis in the tongue site, the two methods 
disagreed. Importantly, differences between the qPCR 
assay and the high-throughput sequence data methods 
may have been the cause of discrepancies. The qPCR assay 
depends on primers that may not be entirely specific to the 
intended organism. The high-throughput sequence data, on 
the other hand, is affected by the reference dataset used to 
classify the sequences. Misclassifications or OTUs not hav- 
ing the required taxonomic resolution would cause discrep- 
ancies when compared with qPCR data. 

Marked differences in a-diversity between health and 
disease consortia were not detected (Figure 1). These re- 
sults are not in agreement with findings reported in two 
recent high-throughput sequencing studies that re- 
ported significant elevations in community a-diversity 
under health and advanced periodontal disease [25,27]. 
It is likely that discrepancies in the sampling technique 
(i.e., curette versus paper-point methods), sampling location, 
and/or methods of DNA extraction can explain differences 
in a-diversity measures between this research and previous 
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studies. Other issues that might have also influenced 
the conclusions could be related to disease activity in 
either the study reported here or others. In this regard, 
even if periodontal pockets are identified, this does not 
mean a priori that the disease is in an active state. In such 
cases, it is also possible that the microbial communities 
found in patients with active versus inactive periodontitis 
might actually be different from one another. Hence, 
further and more detailed study is probably needed to 
consider active versus inactive disease when taking 
microbiological samples. Indeed, it might even be possible 
that a particular shift in the microbiome could be used as 
a predictive marker for the development of active disease 
or for the presence of inactive disease. For a better un- 
derstanding of a-diversity differences, future studies 
should also consider sampling both periodontal pockets 
and healthy tissue from the same patient, controlling for 
the between-patient variability observed here. 

Conclusions 

For the first time, we demonstrate that the tongue dor- 
sum and SupG sites can be used as alternative sampling 
sites for detection and enumeration of bacterial bio- 
markers associated with CP. Use of plaque from these 
sites, as opposed to sampling periodontal pockets has 
important clinical implications. The tongue is an easily 
accessible body site and sampling the tongue dorsum 
for plaque collection is convenient and less invasive. 
This means that others with an interest in examining, for 
example, correlations between CP and certain non-oral 
disease (for example, respiratory disease) could conceivably 
perform meaningful microbiological analyses even if 
they have no training in the dental clinical sciences. It 
is also noteworthy that even within dental research settings, 
ongoing or post-treatment assessments of the micro- 
bial population could be carried out by almost any 
technologist thus reducing the costs of investigation 
since dentists or dental hygienists would not necessar- 
ily be needed for the procedure. Thus, the spectrum of 
those studying oral health and general health correla- 
tions and interactions could be widened considerably, 
at least where microbiological parameters are needed. 
Currently, 27 CP indicator organisms identified in this 
study are being quantified in periodontal and case-matched 
control samples via qPCR to identify the strongest disease 
biomarkers for both SupG and tongue sites. 

The findings reported here add to knowledge of oral 
microbial communities associated with periodontitis. 
Of classic and novel periodontal pathogens identified 
to date, the role of F. alocis as a diagnostic marker for 
CP is noteworthy. Based on results from this study, as well 
as data presented by others, we propose that F. alocis be 
included as a novel member in Socransky s Red complex of 
periodontal pathogens [22]. Because chronic periodontal 



destruction was attributed predominantly to proliferation 
of Gram-negative pathogens, inclusion of F. alocis in the 
Red cluster also changes this classic perspective. 

Methods 

For this cross-sectional study, 24 patients were recruited 
through the University of Toronto Graduate Periodontology 
and Undergraduate Clinics. The study was approved by the 
University of Toronto Research Ethics Board (REB 23872), 
and informed consent to participate in the study was ob- 
tained from all patients. Eleven healthy and 13 patients 
with advanced CP (confirmed by clinical and radiographic 
examinations) consented to participate in this investigation 
(Table 1). Patients with systemic disease and/or those who 
had used antibiotics or had undergone a professional dental 
cleaning within the last six months were excluded from the 
study. Both SupG and SubG plaque samples were collected 
from interproximal sites between teeth 16/15 and 15/14 as 
well as 26/25 and 25/24 using separate metal curettes. After 
the SupG plaque was collected, the tooth was wiped with 
sterile gauze and the subgingival environment was accessed 
through the periodontal pocket. Samples from each site 
were placed in phosphate-buffered saline (PBS) and then 
subsequently pooled into SubG and SupG for each patient; 
samples were not pooled across patients. Plaque from the 
tongue was collected using a curette by making two midline 
scrapes of the dorsal surface. Plaque was then frozen 
at -80°C until DNA analysis could be carried out. All plaque 
samples used in this study were collected by the same clin- 
ician to limit the variability during plaque sampling. DNA 
was extracted using the BacReady DNA extraction kit util- 
izing a single enzyme lysis system as recommended by the 
supplier (GeneScript USA Inc., Piscataway, NJ, USA). 
Briefly, 1 uL of the plaque mixture was added to 20 uL of 
the Buffer BR-A (Genscript, NJ), mixed thoroughly by pip- 
etting and incubated for 10 min at room temperature. 
These reactions containing template DNA from lysed cells 
were used for PCR as described below. To validate the use 
of BacReady DNA kit for DNA extraction, we employed a 
mock bacterial community that comprised equal numbers 
of cells of five Gram-negative and Gram-positive species that 
are found in the oral cavity. These included Actinomycetes 
actinomycetemcomitans, P. gingivalis, T. denticola, T. 
forsythia, and Streptococcus mutans. Lysis mixtures 
were used for PCR amplification using species-specific 
primers (not shown), whereas genomic DNA obtained 
from pure cultures was used as positive controls. Using 
agarose gel-electrophoresis, we obtained the expected 
amplicons for all five species (data not shown) using 
the BacReady kit for cell lysis, which was used to justify 
our DNA extraction protocol. 

A 16S rRNA gene library of approximately 200-bp V3 
16S rRNA genes was generated using modified bacterial 
341 F and 528R primers [43]. Three PCR amplifications 
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were performed for each sample using 50- (il reaction 
volumes in order to minimize potential amplification 
bias [44]. Each reaction mixture contained 25 pmol of 
each primer, 200 \iM of each dNTP, 1.5 nM MgCl 2 and 
1 U of Phusion Taq polymerase (New England Biolabs, 
Ipswich, MA, USA), and 10 ng of template DNA. Am- 
plifications were carried out as follows: denaturation at 
95°C for 5 min, 20 cycles at 95°C for 1 min, 50°C for 1 min, 
and 72°C for 1 min, followed by an extension step at 72°C 
for 7 min. A 2% agarose gel was used to separate products 
from primers; the correct band was recovered using a QIA- 
quick Gel Extraction Kit as recommended by the manufac- 
turer (Qiagen, Valencia, CA, USA). After measuring sample 
concentration using a NanoPhotometer (Implen, Denmark), 
10 ng/(il aliquots of PCR amplicons were combined prior to 
paired-end sequencing with individually indexed samples 
[43]. Sequencing was performed on an Illumina MiSeq 
(2 x 151 bp using a 300 cycle reagent kit). All sequences 
were deposited into the European Nucleotide Archive 
under the project accession number PRJEB6047. 

Sequence data were analyzed in four parallel AXIOME 
(vl.6) analyses (i.e., SubG samples, SupG samples, tongue 
samples, and all three sites combined). All analyses were 
conducted using the same configurations. Paired-end reads 
were assembled using PANDAseq [45], then assembled 
sequences were clustered into operational taxonomic 
units (OTUs) at 97% identity using cd-hit-est (v4.5.4) 
[46]. Multiple sequence alignment of representative se- 
quences was completed with PyNAST (vl.2) via QIIME 
(vl.6) and FastTreeMP (v2.1.3), with default parameters, 
to generate a phylogenetic tree of OTUs to calculate 
UniFrac distances [47-50]. Classification was completed 
using RDP (v2.2) via QIIME with a confidence cutoff of 
0.8 [48,51]. Classifications were rank- flexible, meaning 
each OTU was only classified down to the lowest taxo- 
nomic rank that had 80% or greater posterior probability. 
A merged Greengenes (October 2012 revision) [52,53] 
and OSU CORE oral database (9 February 2012 revision) 
[54] was used for classification. The final OTU table was 
generated by QIIME. UniFrac-based principal coordinates 
analysis (PCoA) which was performed through QIIME. 
Indicator species analysis [29], Bray-Curtis PCoA ordina- 
tions, and multiresponse permutation procedure (MRPP) 
analyses [55], based on nonmetric multidimensional scaling 
(NMDS) plots using Bray-Curtis distances, were conducted 
by the AXIOME pipeline [40]. Alpha-diversity plots were 
generated by Phyloseq (vl.3.14) [56]. 

Validation of CP indicator species using quantitative 
Real-Time PCR (qPCR) 

Forty nanograms of DNA from each plaque sample were 
amplified using qPCR with species-specific primer-pairs for 
Tanerella forsythia (Tf: 5 -GGATTGACCACCGGCGAAG 
ACA-3' and 5 -CGGACACGACGGTTACTCAAATGG-3 '), 



P. gingivalis (Pg: 5 -CTTGACTTCAGTGGCGGCAG-3 ' 
and 5 -AGGGAAGACGGTTTTCACCA-3 '), and Filifactor 
alocis {Fa: 5 ' -TCGGTGCCGAAGTTAACACA-3 ' and 
5 -GGGTCCCCGTCAATTCCTTT-3 ) as well as uni- 
versal primers (5 -TCCTACGGGAGGCAGCAGT-3 ' and 
5 -GG ACTACCAGGGTATCTAATCCTGTT-3 '), to quan- 
tify total bacteria. Each reaction was run in duplicate using 
the Stratagene Mx3000P cycler (Agilent Technologies, 
Santa Clara, CA, USA), and reciprocal values of the mean 
normalized Ct values (Ct indica tor species/ Ct tota i bacteria) were 
obtained. Significant differences between healthy cohorts 
and those with CP within each of the respective locations 
in the oral cavity were calculated based on use of the non- 
parametric Mann- Whitney test with GraphPad Prism 6 
Software (La Jolla, CA, USA) and outliers removed using 
Grubbs Outlier Test (extreme studentized deviate). 

Additional files 



Additional file 1: Sample and Patient Metadata. This includes sample 
ID, age, weight, health status, and amount patient smokes and drinks. 

Additional file 2: UniFrac Principal Coordinates Analysis (PCoA) 
ordination. PCoA ordination of UniFrac distances from a phylogenetic 
tree created by FastTree. A: All samples (SubG, SupG, and Tongue). 
B: SubG samples. C: SupG samples. D: Tongue samples. 

Additional file 3: Multi-Response Permutation Procedure (MRPP) 
results of Health and Disease. MRPP results for each site using 
Bray-Curtis dissimilarity measures. The A statistic is a measure of 
chance-corrected within-group agreement. The closer this value is to 1, 
the stronger the agreement within a group. The J statistic is a 
sample-size dependent measure of the separation between the health 
and disease groups, with more negative J meaning a stronger separation 
between groups. 

Additional file 4: Relative abundance of genera between healthy 
and disease samples for separate SubG, SupG, and Tongue samples. 

The top ten most abundant genus level classifications with significant 
differences between health and disease are shown. Significance was 
assessed with a Mann-Whitney test (P <0.05). Disease abundances are 
shown in pale red, and health abundances in green. "Unci." indicates the 
group was not classified down to genus, and the lowest level classification 
available is given. A: SubG samples. B: SupG samples. C: Tongue samples. 

Additional file 5: Relative abundance of species between healthy 
and disease samples. The top ten most abundant species level 
classifications with significant differences between health and disease are 
shown. Significance was assessed with a Mann-Whitney test (P < 0.05). 
Disease abundances are shown in pale red, and health abundances in green. 
"Unci." indicates the group was not classified down to species, and the lowest 
level classification available is given. A: All samples (SubG, SupG, and Tongue). 
B: SubG samples. C: SupG samples. D: Tongue samples. 

Additional file 6: List of species significantly increased under health 

or disease. Spreadsheet file containing significance results comparing 
relative abundance of taxa under health and disease for SubG, SupG, and 
Tongue sites. Significance was assessed with a Mann-Whitney test (P <0.05). 

Additional file 7: Dufrene-Legendre indicator species analysis. 

Spreadsheet file containing the detailed indicator OTU lists for the 
combined, SubG, SupG, and Tongue analyses. All indicator OTUs with 
P< 0.05 are included. 
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